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Abstract 

The probability that a generic real, complex or quaternionic two-qubit state is separable can 
be considered to be the sum of three contributions. One is from those states that are absolutely 
separable, that is those (which can not be entangled by unitary transformations) for which the 
maximal concurrence over spectral orbits {Cmax) is zero. The other two contributions are from the 



states for which Cmax G (0, and for which Cmax G [5) !]• We have previously (arXiv:0805.0267) 
found exact formulas for the absolutely separable contributions in terms of the Hilbert-Schmidt 
metric over the quantum states, and here advance hypotheses as to the exact contributions for 
Cmax G [5,1]. A crucial element in understanding the two contributions for Cmax > is the 
nature of the ratio (R) of the Cmax-parameterized separability function for the complex states 
to the square of the comparable function for the real states-both such functions having clearly 
displayed jump discontinuities at Cmax = \- For Cmax £ (O)^]) the ratio R appears to be of 
the form 1 + kCmax, except near Cmax = \i while for Cmax £ there is strong numerical 

evidence that it equals 2 (thus, according to the Dyson- index pattern of random matrix theory). 
Related phenomena also occur for the minimally-degenerate two-qubit states and the qubit-qutrit 
states. Our results have immediate application to the computation of separability probabilities 
in terms of other metrics, such as the Bures (minimal monotone) metric. The paper begins with 
continuous embeddings of the separability probability question in terms of four metrics of interest, 
using " generalized Peres-Horodecki conditions" . 
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I. INTRODUCTION 

One possibly productive strategy to pursue when confronted with an apparently in- 
tractable problem, is to embed it in some broader class of problems. Doing so, hopefully, 
may lead to new insights and progress, including ones regarding the original (smaller) prob- 



lem. In the first of the two basic parts of this paper (sees. II VII), we adopt such embedding 
strategies for the task of determining the probabilities-with regard to a number of metrics 
of quantum-mechanical interest-that certain generic forms of 2 x 2 or 2 x 3 quantum sys- 
tem are separable [Ill2l|3llH|5l|6|,[7|,|8]. (As computers presumably grow more powerful, 
these readily-formulated, but high- dimensional [9, 15, ...] and high-degree [e. g., quartic] 
problems may eventually lose their apparent present-day intractability (cf. [9l[T0])-much as 
did the famous four-color planar map theorem of Appel and Haken [HI [12]. Nevertheless, it 
would certainly be appealing to address these problems with more theoretical understanding 
than is required by "brute force" computation (cf. [131 IH [151 US])-) 
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In the second basic part (sees. VIII X), building upon our recent work in [T7], we attempt 
to gain insight-using manifest relations to random matrix theory-into the very same sepa- 
rability probability questions by determining the nature of certain eigenvalue-parameterized 
separability functions. These are expressed as univariate functions of the maximal concur- 
rence over spectral orbits. 

The fundamental question being addressed here of determining the probability that a 
generic bipartite quantum state is separable or not was first raised by Zyczkowski, Horodecki, 
Sanpera and Lewenstein in a pioneering, much-cited 1998 paper [T]. As motivation they 
wrote: "One of the fundamental questions concerning these subjects is to estimate how 
many entangled (disentangled) states exist among all quantum states. More precisely, one 
can consider the problem of quantum separability or inseparability from a measurement 
theoretical point of view, and ask about relative volumes of both sets. There are three main 
reasons of importance in this problem. The first reason, of some philosophical implication, 
may be contained in the questions Is the world more classical or more quantum? Does it 
contain more quantum-correlated (entangled) states than classically correlated ones? The 
second reason has a more practical origin. Analyzing some features of entanglement, one 
often has to rely on numerical simulations. It is then important to know to what extent 
entangled quantum states may be considered as typical. Finally, the third reason has a 
physical origin. The physical meaning of separability has recently been associated with the 
possibility of partial time reversal" [H p. 883]. 

In sec. |Tl| in the first basic part of the paper, we analyze the cases of generic 9-dimensional 
real and 15-dimensional complex two-qubit systems. For our calculations, we utilize the 
Euler-angle parameterizations of the real (developed by S. Cacciatori [HI App. A]) and of 
the complex 4x4 density matrices (p) p^, as well as the Tezuka-Faure (TF) procedure 
[T9t 120] for generating low- discrepancy sets of (9- and 15-dimensional) points. These points 
are employed for quasi-Monte Carlo numerical integration with respect to the product of 
the (6- or 12-dimensional) Haar measure over the Euler angles and (3-dimensional) metric- 
specific measures over the eigenvalues of the density matrices. In sec. |III[ we turn our 
attention to parallel "continuous embedding" analyses pertaining to the 14-dimensional 
(rank-3) boundary of the 15-dimensional generic complex 4x4 density matrices. 



In sec. IV, we make use of the S'f/(6)-based Euler-angle parameterization of the 35- 
dimensional generic complex quhit- qutrit 6x6 density matrices [21] sec. XI] to investigate 
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the corresponding rank-6 and rank-5 problems. In sec. |V} we investigate formally extending 
the range of our basic parameter (Q;)-used in forming convex combinations-from beyond [0,1] 



to [—00,00]. In sec. VI, we depart from our initial paradigm, employing the parameter a, 
and evaluate the separability probabilities of the generic complex and real two-qubit states 
for which the entanglement measure concurrence (C) [221 ES] is less than some threshold. 
(There we observe some interesting behavior, involving the intersection of the curves for 
different metrics. For C = 1, we obtain the usual separability probabilities.) 



The concept of maximal concurrence (Cmax) over spectral orbits 

Cmax = max(0, Ai - A3 - 2a/Vm), (Ai > A2 > A3 > A4) (1) 

(a quantity which can not be increased under unitary transformations) of a two-qubit density 
matrix (p), where the A's are the ordered eigenvalues of p, is used in the second basic set of 



analyses of the paper (sec. VIII). There, we importantly add to certain findings [17] concern- 
ing the strong goodness-of-fit to two-qubit eigenvalue-parameterized separability functions 
(ESFs) of piecewise functions of Cmax- We find evidence of adherence over a half-domain 
Cmax £ [|, 1] to a Dyson-index pattern both for the generic rank-4 (first investigated in 
|17j ) and generic rank-3 (as found here [sec. VIII E 1| ) real and complex two-qubit states. 



We, further, undertake an analogous examination of: (a) the generic rank-5 qubit-qutrit 
states in sec. VIII F, observing interesting jump discontinuities in the ESFs; and (b) the 



generic full rank qubit-qutrit states in sec. |VIII G[ where, again, the Dyson-index pattern 
appears to emerge over a restricted domain Cmax & [IjI]- Additionally, in sec. 



IX A 



we 



are able to present new simple exact results pertaining to certain components of the de- 
sired Hilbert-Schmidt separability probabilities. In these regards, let us draw the reader's 



attention, particularly, to (the titular) Figs. 34 and 35 



Since when we had earlier addressed the issue of two-qubit separability probabilities in 
terms of diagonal- entry-paiameterized separability functions (DESFs) [7], we found appar- 
ently total agreement with Dyson-index behavior, the need (remaining unmet) to reconcile 
these two forms of Dyson-index patterns (full and partial) is obvious. 

Remarks relevant to the two primary sets of analyses-which share the use of concurrence 



and are devoted to the determination of separability probabilities-are given in sees. VII and 

Kl 
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A. Separability functions 



Let us state here that the concept of a separability function-hoth in its eigenvalue- 
parameterized (ESF) and diagonal-entry-parameterized (DESF) forms-has been developed 
in order to reduce the intrinsically high dimensionalities of the generic separability proba- 
bility questions. By integrating over the majority of parameters-for example, Euler angles 
or off-diagonal entries-one reduces-the problems (at least, in the two-qubit case) to ones 
of (only!) a three-dimensional nature. It also appears possible to further reduce the three- 
dimensional problems to single-dimensional ones by finding an appropriate parameter-known 
to be the ratio of the product of the 11- and 44-diagonal entries to the product of the 22- 
and 33-diagonal entries in the DESF-case, and apparently (as numerics strongly indicate) 
the maximal concurrence over spectral orbits in the ESF-case. 

B. Metrics employed 

The metrics of quantum-mechanical interest that we utilize to form [via their Rie- 
mannian volume elements) measures over the quantum systems are the (Euclidean or 
flat, non-monotone [23]) Hilbert-Schmidt [I3j, and three monotone metrics [2n]^the (min- 
imal monotone) Bures [14], Wigner-Yanase [27] and Kubo-Mori [28] metrics. (We also 
attempted to include the information-theoretically significant monotone "quasi-Bures" 
[Grosse-Krattenthaler-Slater] metric [5l|29l|30], which yields the minimax/maximin asymp- 
totic redundancy for universal quantum coding, but encountered some initial, at least, nu- 
merical difficulties in this regard.) 

C. Prior conjectures 

In [7], we were led by a combination of numerical and theoretical [Dyson-index-related] 
arguments, involving diagonal- entry-paiameteiized separability functions (DESFs), to con- 
jecture that the Hilbert-Schmidt separability probabilities are, respectively, ~ 0.242424 
for the generic complex two-qubit systems and ^ ~ 0.470588 for its real counterpart. (The 
supporting evidence appeared particularly strong for the ^ figure.) Also, it has been fur- 
ther conjectured that for the generic quaternionic two-qubit systems, the corresponding 
probability is ^ 0.0733389 [SD eq. (15)] [SI p. 25]). 
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We had also earlier advanced in ^ Table VI], the "silver mean" (that is, — 1) con- 
jectures that the generic complex two-qubit Bures separability probability is 

^ 1680(v^- 1) _ Q Q7333gg /2) 

Sep/ Bures ^8 ^ ' 

and the corresponding Kubo-Mori analogue is 

= '^^^2^"'^ - 0.035398. (3) 

(Additionally, for the "average monotone metric"-not employed in this paper-it was con- 
jectured in [1] that the associated separability probability is ^^^^^^"^'^^ ^ 0.0475329. Fur- 
ther still, the Wigner-Yanase separability probabihty was hypothesized to equal the ratio 
of '^^^~^'> to the not-yet-determined Wigner-Yanase volume of the generic [entangled and 
separable] complex two-qubit states.) 

In [71 sec. X], again studying the corresponding DESFs, the conjectures were put forth 
that the generic real and complex quhit-qutrit Hilbert-Schmidt separability probabilities are, 
^ ^ 0.150235 and (agreeing very closely with the numerics) ~ 0.0266889, respectively. 



II. GENERIC FULL-RANK REAL AND COMPLEX TWO-QUBIT CASES 
A. First set of constraints convex combinations of determinants of p and ppT 

1. Real two-qubit density matrices 

In Fig. [l]we show as a function of a G [0, 1] the probabilities {Pmetrid^^)) y terms of the 
four metrics under consideration, that for a generic (9-dimensional) real two-qubit system 

a\ppT\ + {I- a)\p\>0. (4) 

Here ppt is the partial transpose of p and |ppt|, its determinant. Of course, here and 
throughout we incorporate into our analyses, the original, notable Peres-Horodecki necessary 
and sufficient conditions for separability in terms of the nonnegativity of ppt The 
partial transpose of a 4 x 4 density matrix can have at most one negative eigenvalue, so the 
condition \ppt\ < is fully equivalent to ppt having a single negative eigenvalue |34j. (Also, 
obviously, the nonnegativity condition |p| > is always satisfied.) 
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FIG. 1: Curves generated by enforcement of determinant-based constraint Q for the generic 9- 
dimensional real density matrices. The order of dominance of the four metric-specific curves is 
given in 17,100,000 Tezuka-Faure 9-dimensional points were employed in the quasi-Monte 

Carlo numerical integration. The values at a = 1 are the embedded (conventional) separability 
probabilities. 

The order of dominance of the four monotonically-decreasing curves in Fig. [T| as well as 
all the other analogous curves below, turns out-with the important exception of those in 
sec. 



VI B, where we observe intersecting behavior-to be 



Hilbert-Schmidt > Bures > Wigner-Yanase > Kubo-Mori. 



(5) 



This, of course, will imply that the associated two-qubit separability probabilities (corre- 
sponding to a = 1) adhere to the same ordering. Since the Bures metric is also the minimal 
monotone metric, it is not surprising that it is extremal among the three monotone met- 
rics under consideration. In estimating these curves, as well as all others displayed below 



involving a-except Fig. 14 -we subdivided the unit interval a G [0,1] into one thousand 
subintervals. 

We can fit the Hilbert-Schmidt curve in Fig. [l] rather well-the integral over a E [0, 1] 
of the sum of squares of the differences being only 0.00026557-while exactly achieving the 
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Of 



FIG. 2: Curves generated by enforcement of determinant-based constraint Q for the generic 15- 
dimensional complex density matrices. The order of dominance of the four metric-specific curves 
is given in ([s]). 15,400,000 Tezuka-Faure 15-dimensional points were employed in the numerical 
integration. 



conjectured separability probability of ^, with the simple function 



HS ' 



a) 
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8 + 9^/a 



(6) 



2. Complex two-qubit density matrices 



In Fig. M we analogously show as a function of a G [0,1], the four probabilities 



{Pmetric^ i'^)) ^ generic (15-dimensional) complex two-qubit system that the inequality 
constraint Q is satisfied. We can fit the Hilbert-Schmidt curve here very well-the integral 
over a G [0, 1] of the sum of squares of the differences being only 0.00030092-while achieving 
our conjectured separability probability of ^ [7] with the function 



-^complex , 
HS ' 



a) 



(7) 



where c = 8 + 2-\/66. Also, the square of the real counterpart ^ does provide a close 
fit to the Hilbert-Schmidt curve in Fig. [2] (However, our conjectured complex two-qubit 
separability probability of 



^ ^ 0.242424 is not equal to the square, J§ ^ 0.221453, of the 
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conjectured real separability probability of so conformity to a Dyson-index pattern is 
not total.) 

Additionally, we can very well fit the Bures curve in Fig. |2] and our corresponding con- 
jectured "silver mean" separability probability (|2]) |1] by the function (of the seventh root 
of a) 

-ocomplexf \ 1680(1 — \/2) , , 

^Bures l«J " ^ T^' l^J 



(^4y/l05 (-1 + V2) (^ - 1) + 7r4^^ 



(The sum-of-squares measure of fit between the two curves is 0.000739208. However, the 
[exact] square root of ([sj-deviating substantially from a Dyson-index-like pattern-does not 
at all provide a close fit [as in the HS case] to the real Bures counterpart in Fig. [l|) 



B. Second set of constraints convex combinations of minimum eigenvalues of p 
and PPT 

1. Real two-qubit density matrices 

In Fig. [3] we show as a function of a G [0, 1], the probabilities for a generic (9-dimensional) 
real two-qubit system that 

O^^PTmin + (1 - 0!)Xmin > 0, (9) 

where the subscript min denotes the smallest of the corresponding four eigenvalues (A) of 
either p or ppT- (As noted, having all eigenvalues nonnegative is fully equivalent to having 
a nonnegative determinant for the partial transpose of a 4 x 4 density matrix |3l]. The 
entanglement measure negativity is equal to max[0, — 2Ap7;^;^] [HI p. 401]. The Hilbert- 
Schmidt distance of an entangled state to the set of all partially transposed sets can be 
expressed as a function of the negative eigenvalues of the partial transpose of the entangled 
state [35].) 



2. Complex two-qubit density matrices 

In Fig.|4]we show as a function of a G [0, 1], the probabilities for a generic (15-dimensional) 
complex two-qubit system that the inequality constraint ^ holds. One can fit the complex 
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FIG. 3: Curves generated by enforcement of minimum-eigenvalue-based constraint ^ for the 
generic 9-dimensional real density matrices. The order of dominance of curves is given in 
32,700,000 Tezuka-Faure 9-dimensional points were employed in the numerical integration. 




I I I I I I I I I I I I I I I I I I I I ^ 

0.2 i)A CIj6 O.B 1.0 

FIG. 4: Curves generated by enforcement of minimum-eigenvalue-based constraint ^ for the 
generic 15-dimensional complex density matrices. The order of dominance of curves is given in 
18,450,000 Tezuka-Faure 15-dimensional points were employed. 
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HS curve rather closely by the square of the corresponding real HS curve, particularly so if 
one adds a small linearly increasing correction of the form ^ to this square. 



C. Third set of constraints determinants of convex combinations of p and ppx 

1. Real two-qubit density matrices 

In Fig.jsjwe show as a function of a G [0, 1], the probabilities for a generic (9-dimensional) 
real two-qubit system that the positive "twofold partial" transpose condition, 

\appT + (1 - a)p\ > 0, (10) 

holds. The Hilbert-Schmidt curve is highly linear in character. The line 1 — |y closely 
approximates it, as well as reproducing the conjectured separability probability of ^. 



2. Complex two-qubit density matrices 

In Fig. [g] we show as a function of a G [0, 1], the four probabilities {Pm^ric^i^)) ^ 
generic complex two-qubit system that the inequality ^ holds. 

In computing Figs. |5] and [6| we solve the quartic equation \appT + (1 — a)p\ = and 
assume that there can not be more than one solution a G [0,1]. (Numerically, this did 
appear to be the case, except for some isolated instances in which two positive [essentially 
identical] roots both very close to zero were found.) 

Another possible generalized Peres-Horodecki condition-that is, that the minimum eigen- 
value of appT + (1 — a)p be nonnegative-appeared to be considerably more problematical 



(time-consuming) than ( 10 ) to investigate, though we do, in fact, implement such a condition 



for the generic complex qubit-qutrit systems (sec. IV) and in generating Fig. 14 

Let us note that all the curves displayed so far in this communication appear to correspond 
to convex functions, but for the last two Hilbert-Schmidt curves. 



D. Comparison of metric-specific curves for the first three sets of constraints 

In Fig. [Tjwe show in a single plot, the three Hilbert-Schmidt curves plotted above (one 
per figure) for the generic real 9-dimensional two-qubit systems, while in Fig. [sjwe show in 
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FIG. 7: The curves corresponding to the Hilbert-Schmidt metric plotted above for the three differ- 
ent constraints for the generic real 9-dimensional two-qubit systems. The quasi-linear curve based 
on constraint (10) dominates that based on ([9|, which dominates that based on Q. 



a single plot, the three Bures curves plotted above for the generic complex 15-dimensional 
two-qubit systems. In both these figures the order of dominance of the curves is the same-the 
curve based on the constraint (10) dominates that based on ([9]), which, in turn, dominates 
that based on (|4]). 

Of course, we find in both of these figures that the three curves have common points-of- 
intersection at a = 1, corresponding to the "ordinary" separability probability (as well as 
a = 0). 



III. GENERIC RANK-3 COMPLEX TWO-QUBIT CASE 



The "twofold" volume-to-area-ratio theorem of Szarek, Bengtsson and Zyczkowski [36] 
allows us to immediately extend our conjecture [7] of Hilbert-Schmidt separability proba- 
bility of for the 15-dimensional generic complex two-qubit states to the fully equivalent 
conjecture that the HS separability probability of the states on the 14-dimensional boundary 
is one-half oi this, that is, In Figs, [o] and 10 we show our corresponding estimation of 
the a-separabilities based on certain obvious modifications of the determinant constraint ^ 
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FIG. 8: The curves corresponding to the Bures metric plotted above for the three different con- 
straints for the generic complex 15-dimensional two-qubit systems. The curve based on constraint 
(10) dominates that based on ([9|, which dominates that based on Q. 



and the minimum eigenvalue constraint ([o]). That is, rather than using the (zero) determi- 
nant of the rank-three density matrix, we use its generically nonzero 3x3 principal minor. 
Further, rather than using the minimum (zero) eigenvalue, we employ the minimum of the 



generically three nonnegative eigenvalues. In Fig. [IT] we display the rank-3 a-separability 
probability estimates based on the application of the constraint (10). Here we notice some 
unusual behavior near a = due to the degeneracy (zero determinant) of a rank-3 two-qubit 
(4 X 4) density matrix. 



IV. GENERIC FULL-RANK COMPLEX QUBIT-QUTRIT CASE 

Of the three distinct sets of constraints considered in the two-qubit case, only ^ seemed 
immediately adoptable to the qubit-qutrit case associated with 6x6 density matrices. In 
our computations, we now employ the associated SU{Q) Euler-angle parameterization [2n 



sec XI]. In Fig. 12 we show the corresponding plot. 

Further, by specifically checking nonnegativity at each value of a = ... 1, we were 
able to enforce the condition that the minimum eigenvalue of the matrix convex combination 
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FIG. 9: Curves generated by enforcement of modified (3x3 principal minor) constraint Q for 
the generic 14-dimensional rank-3 complex density matrices. The order of dominance of curves is 
given m 7,300,000 TF 14-dimensional points were employed. 




FIG. 10: Curves generated by enforcement of modified (minimum-eigenvalue-based) constraint ([9]) 
for the generic 14-dimensional rank-3 complex density matrices. The order of dominance of curves 
IS given m 10,300,000 TF 14-dimensional points were employed. 
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FIG. 11: Curves generated by enforcement of constraint (10) for the generic 14-dimensional rank- 



3 complex density matrices. The order of dominance of curves is given in 1^. 4,300,000 TF 
14-dimensional points were employed. 

apPT + (1 — «)p be nonnegative. (Nonnegativity of the determinant of the partial transpose 
is no longer equivalent-as it is in the two-qubit case [3l]-to having no negative eigenvalues, 
since two negative eigenvalues yields a positive determinant.) The corresponding plot is 



displayed in Fig. 13 



We can "sandwich" the Hilbert-Schmidt curve in Fig. 12 between two curves, correspond- 



ing to exact squares, both of which yield the conjectured separability probability of jj|g. 



These functions are 



64 



((-8 + V2398) + 



and 



64 



((-8 + V2398)a + 8) 



2 • 



(11) 



V. EXTENDING RANGE OF a-PARAMETER 



We have, so far, considered our primary variable a as extending over the unit interval 
[0,1]. However, it appears quite interesting and possibly more natural to formally view its 



range as the real line [— cx3,cx3]. In a further analysis, we developed a plot (Fig. 14) over 
a G [— |, ^], of the estimated a-probability that the 4x4 matrix appT + {1 — Ci)p, where p 
is a generic complex two-qubit density matrix, has all its four eigenvalues nonnegative. 
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FIG. 12: Curves generated by enforcement of (minimum-eigenvalue-based) constraint ([9]) for the 
generic 35-dimensional complex 6x6 density matrices. The order of dominance of curves is given 
in 30,650,000 TF 35-dimensional points were employed. 




FIG. 13: Curves generated by verifying the nonnegativity of the minimum eigenvalue of appT + 
(1 — a)p for the generic complex qubit-qutrit systems. The order of dominance of curves is given 
in 1,300,000 TF 35-dimensional points were employed. 
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T : 1 



-2 



FIG. 14: Curves generated by verifying the nonnegativity of the minimum eigenvalue of appT + 
(1 — a)p for the generic complex two-qubit systems, using the extended interval a E [— fj The 
order of dominance of curves is given in ([5|. 1,650,000 TF 15-dimensional points were employed. 

VI. CONCURRENCE-RELATED ANALYSES 
A. Generalized Peres-Horodecki conditions 

In all the analyses reported above, the nonnegativity convex combination constraints 
("generalized Peres-Horodecki conditions") utilized, have been expressed either in terms of 
the determinant or the minimum eigenvalue of Ppt- In the two-qubit case, we have also 
been able to investigate similiarly-motivated conditions using, in conjunction, the maximal 
concurrence over spectral orbits (jTl) [24J of a two-qubit density matrix (p), and its concurrence 



(Here, the A's are the ordered eigenvalues of p and the r^'s are the ordered eigenvalues 



denotes conjugation. Throughout the reminder of the paper, the symbol cr-consistently 
with our previous notation- will denote a "separability function, and not a Pauli matrix.) 




C = max(0, T]i-ri2-r]3- rj^) 



iVl >V2>V3> Vi)- 



(12) 
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FIG. 15: Implementation of concurrence-based constraint (13) for the generic 9-dimensional real 



two-qubit states. The order of dominance of curves is given in ([5]). 9,250,000 TF-points were 
employed. 

The corresponding constraint we employ is 



- aC + {1 - a)Cmax > 0. 



(13) 



Since Cmax > C, the constraint holds trivially for a G [0, |]. In Figs. 15 and 16, we show 
for the /iail/-interval a G [|, 1], the curves for the generic real and complex two-qubit states, 



respectively, based on while in Figs. 17 and 18, we display the corresponding plots for 



the generic rank-3 real and complex two-qubit states, respectively. 



B. Separability probabilities as functions of concurrence— intersecting curves 

In this section, we depart from the basic paradigm so far employed in first basic part 
("generalized Peres- Horodecki conditions") of the paper, in which we use convex combina- 
tions of quantum-mechanical terms to form nonnegativity constraints. 

Now, we simply estimate-again, with respect to the four metrics in question-the separa- 
bility probability of two-qubit states for which the concurrence C is less than some threshold 



Cq. We show our results in Figs. 19 and 20 The generic rank-3 counterparts of these two 
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FIG. 16: Implementation of concurrence-based constraint (13) for the generic 15-dimensional com- 



plex two-qubit states. The order of dominance of curves is given in (H. 6,100,000 TF-points were 
employed. 
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FIG. 19: Separability probabilities of generic complex two-qubit states having concurrence less than 
or equal to C. The Hilbert-Schmidt (blue) curve intersects the other three. 21,100,000 TF-points 
were employed. 
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FIG. 20: Separability probabilities of generic real two-qubit states having concurrence less than 
or equal to C. The Hilbert-Schmidt (blue) curve intersects the other three from below, all near 
a = 0.12. There were 7,800,000 TF-points employed. 



figures are given in Figs. [21] and [22} In all four of these cases, the Hilbert-Schmidt curve 
intersects the curves for the three monotone metrics from below. In this regard, it has been 
noted by Bengtsson and Zyczkowski that the "Bures measure is concentrated at the states 
of higher [than the Hilbert-Schmidt] purity" [151 p. 356], since {Ttp^)hs < {^"^P^) Bures- Our 
(intersecting) results in this set of concurrence-based analyses is clearly consistent-but now 
taking a separability-related form-with that assertion. 



VII. REMARKS 



Our motivation in undertaking the first principal part of this study reported above has 
been to examine whether it might be feasible to shift the question of determining the two- 
qubit separability probabilities with respect to various metrics of quantum-mechanical in- 
terest to the (perhaps more tangible, addressable) question of characterizing the curves that 
interpolate between such separability probabilities and the (unit) probabilities that a two- 
qubit state is either separable or entangled (cf. [13], E]). We intend to study the curves 
generated in still greater detail, as additional computations render them more precise. (The 
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FIG. 21: Separability probabilities of generic rank-3 complex two-qubit states having concurrence 
less than or equal to C. The Hilbert-Schmidt (blue) curve intersects the other three. 14,000,000 
TF-points were employed. 




FIG. 22: Separability probabilities of generic rank-3 real two-qubit states having concurrence less 
than or equal to C. The Hilbert-Schmidt (blue) curve intersects the other three. 9,750,000 TF- 
points were employed. 
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Tezuka-Faure procedure is not amenable to use of statistical tests, though variants of this 
quasi-Monte Carlo method have been developed that are.) In particular, it would be of 
interest to see if the differences between the curves for the three monotone metrics studied 
could be explained directly in terms of the Chentsov-Morozova functions c(x, y) for those 
metrics |26l[37j. These are (^^^^y and , for the Bures, Wigner-Yanase and 

Kubo-Mori metrics, respectively. (The associated operator monotone functions, /(t), for 
which c{x,y) = yf^^L^ are *+^'^+^ and respectively.) The possible relevance of the 
Dyson-index- ansafe to explain differences between results for the generic real and generic 



complex systems-as in sec. VIII below-should also be examined [3 [T7] . 

Conceiveably, our attempted generalizations here of the Peres-Horodecki conditions and 
introduction of the concept of " a-separability" might prove productive in some manner 
parallel to the well-studied concepts (also based on generalizations/extensions/embeddings) 
of p-Renyi-entropy [38] and of escort distributions [391 



VIII. SEPARABILITIES AS PIECEWISE CONTINUOUS FUNCTIONS OF 
MAXIMAL CONCURRENCE 

A. Objective 

Here, we begin the second basic part of our paper. We importantly amend a certain 
parenthetical remark made in our recent paper [T7|, to the effect that although two-qubit 
dia^'ona/- enir?/-parameterized separability functions (DESFs) had been shown [3, |3T] to 
clearly conform to a pattern dictated by the "Dyson indices" {(3 = 1 [real], 2 [complex], 
4 [quaternionic] ) of random matrix theory, this did not appear to be the case with regard to 
ezg'enfa/we-parameterized separability functions (ESFs). (We remark here that the "value 
of (3 is given by the number of independent degrees of freedom per matrix element and is 
determined by the antiunitary symmetries ... It is a concept that originated in Random 
Matrix Theory and is important for the Cartan classification of symmetric spaces" [HI p. 
480]. The Dyson index corresponds to the "multiplicity of ordinary roots", in the termi- 
nology of symmetric spaces [121 Table 2].) But upon further examination of the extensive 
numerical analyses reported in [17], we found quite convincing evidence that adherence to 
the Dyson-index pattern does also hold for ESFs, at least as regards the upper half-ismge 
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I < Cmax < 1 of the maximal concurrence over spectral orbits ([l]). 

To be specific, it strongly appears that in this upper half-range, the real two-qubit ESF 

3 

is equal to to i^n^^^^H^ and its complex counterpart-in conformity to the Dyson-index 
pattern-proportional to the square of the real ESF, that is, ^ ^^"^""^ . The previously doc- 
umented piecewise continuous ( "semilinear" ) behavior in the lower half-range < Cmax < | 
appeared to lack any particular Dyson-index-related interpretation-which seemed somewhat 
paradoxical in terms of our DESF-findings [7[ [31] . However, we report new insights into this 



problem below (sec. IXC). 



B. Previous ESF findings 

The study [T7] had been devoted to the question of determining for the generic (9- 
dimensional) real and (15-dimensional) complex two-qubit systems, the nature of cer- 
tain trivariate "eigenvalue-parameterized separability functions" (ESFs). These (metric- 
independent) ESFs, it was argued, could substantially assist in the determination of sepa- 
rability probabilities in terms of certain metrics (the Hilbert- Schmidt and Bures being the 
most conspicuous examples). (In [7], DESFs were successfully used in the Hilbert-Schmidt 
case, but they do not seem as useful for the Bures and other montone metrics, the standard 
formulas for which are expressed in terms of eigenvalues, and not diagonal entries.) We 
further investigated in [IT] the possibility that these prima facie trivariate functions of the 
eigenvalues Aj (i = 1, . . . 4) of 4 x 4 density matrices (A4 = 1 — SfAj), were expressible as 
univariate functions 

5f (Ai . . . A4) = cr^^\Cmaxi\i ■ ■ ■ A4)), (14) 

of the maximal concurrence C^ax^given by ([T])-over spectral orbits fl^^ sec. VII] jlHl SI]. 
(At this point in our presentation, let us-motivated by Dyson-index conventions-regard f3 



in (14) only as a notational [dummy variable], not calculational device taking the values 1 



[real], 2 [complex], 4 [quaternionic] . ) 



C. Jump discontinuities 



Our main conclusions in p!7] were that-if the reducibility-to-univariance property (14) 



held, as our extensive numerical evidence appeared to suggest might be the case (being able 
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normalized sep functs. 
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FIG. 23: Joint plot of numerical estimates of ( ^^^^^V and '^"^^"^7^ for Cmax G \h !]• 

to explain almost 99% of the variance [T7t Sec. II.B.l])-the associated real and complex 
univariate functions both had jumps of approximately 50% magnitude at Cmax = |, as 
well as a number of additional discontinuities (remarkably coincident in both the real and 
complex cases) in the lower half-range Cmax G (0, |]. (The joint jumps at Cmax = \ were 
displayed in [T7] in Figs. 2 and 6. We have since found a small programming error that 
caused the two curves in Fig. 2 there to be slightly more misaligned-by ^-than they should 
have been.) Also, both univariate functions appeared to be simply /mear between certain of 
these discontinuities. The upper half-range Cmax ^ [\i l]-in which the univariate functions 
of Cmax took Icsscr values-did not command our attention in [T^ , seeming to be of relatively 
less interest. Our only pertinent observation there was that there did not appear to be any 
discontinuities in that segment. 




D. New Dyson-index-related findings 

Now, in fact, turning our attention more closely to this upper half-range Cmax G 1]; 
we readily find strong evidence for a very interesting Dyson-index-type phenomenon. If we 
normalize our extensive numerical estimates from [17] of <J^{Cmax) and cr^(Cmaz) to both 



equal 1 at the jump discontinuity point Cmax = \^ then a joint plot (Fig. 23) of the latter 
normalized (complex) function versus the square of the former normalized (real) function 
for Cmax € [^,1] remarkably shows no perceptible difference between the two resulting 



curves. (The sample [quasi-Monte Carlo] estimate of cr^(|) is 0.0651586 and that of cr^(|)) 



is 0.1803748.) In Fig. uM we show-on a much finer scale than used in Fig. 23 -the actual 
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FIG. 24: Numerical estimate of 

normalized sep functs. 
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FIG. 25: The two functions in Fig. 23 along with the additional {very closely- fitting) function 

(2 '^CfYiax) • 



(very small) numerically-obtained differences 



a 



a 



2(1: 



(15) 



between them. Of further considerable importance, Fig. [25] is a repetition of Fig. [23| but 
along with the insertion now of the function 



(2 2Cmux) ^(-^ Cxnax) 



(16) 



which we see fits our two estimates very well. Assuming that (16) is the correct form (up 
to the still not exactly-known normalization factor) of o'^iCmax) over Cmax ^ [|) 1]) we can 
estimate the associated contribution to the separability probabilities from density matrices 
corresponding to this half-range to the Hilbert-Schmidt and Bures separability probabilities 
of generic complex two-qubit systems to be 0.0100578 and 0.0194829, respectively. (The 
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real counterparts of these separability probabilities are, then, 0.0254346 and 0.0100578, 



respectively [cf. ([29]), (|30])].) 

Let us further note that our sample estimate of the ratio 



a^il) 0.0651586 



0.18037482 



2.00272 (17) 



is very close (and possibly theoretically exactly equal) to 2. 

Over < Cmax < |, the range of primary interest in [17j, the estimates of the real and 
complex two-qubit separability functions intersect (near Cmax = 0.1812), and appear to have 
linear segments over the same subintervals [T71 Figs. 1, 5, 7]. These features appeared to 
make any immediate application of the Dyson-index pattern problematical in this lower half- 
range. So, the behaviors of the univariate functions cr^'^\Cmax), (/? = 1 [real], 2 [complex]), 
over the two indicated regimes of Cmax seem to be highly distinct. The point Cmax = | 
clearly serves as a point of major behavioral transition, with the lower half-range, then, 
appearing perhaps to be the more theoretically challenging of the two. (We will observe 
what appears to be similarly dichotomous Dyson-index behavior in the qubit-qutrit case 



[Fig. 31 . Perhaps one might view the two regimes as semiclassical and quantum in nature.) 
An outstanding question is what are the specific values of cr^(^) and cr^(^), which we used 



as normalization factors in our analyses above. The nearness to 2 of the ratio (17) may be a 
helpful guide in this regard. In fact, let us take this opportunity to further indicate that in 
our ongoing supplemental analyses-in which we use 5,000, rather than 500 sampling points 



in the interval [0,l]-we have obtained for the ratio (17) the estimate 



^2(1) 0.066663 , , 

- - 1.99605. (18) 



2 0.182752 



This ratio would be exactly 2 if we took for the numerator of (18) the value ^ ~ 0.0666667 
and for the denominator, ^ ^ 0.1825742. We will, in fact, assume these exact values 



in seeking to ascertain in sec. IX, the exact contributions over Cmax £ [-k, 1] to the total 



Hilbert- Schmidt two-qubit generic real and complex separability probabilities. 



E. Rank-three complex and real two-qubit cases 



Here, we straightforwardly apply the same maximal-concurrence ansatz (14) just dis- 



cussed and applied to the full (rank-4) complex and real two-qubit cases, to the minimally 
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iniik-3 



FIG. 26: Eigenvalue-parameterized separability functions, expressed in terms of maximal concur- 
rence, for the minimally degenerate rank-three generic real and complex two-qubit states. The 
complex (red) curve is initially dominant. For each TF-point employed (2,062,400 in the com- 
plex case and 2,331,300 in the real case), separability is checked for 500 equally-spaced values of 



^max ^- There appear to be discontinuities at CI 



ranks 
max 



1 

2- 



degenerate (rank-3) counterparts. The main conceptual point to note is that the formula 
for the maximal concurrence ([T| now degenerates to 



a 



rank— 3 
max 



Ai-A 



3; 



(Ai > As > A3). 



(19) 



In Fig. 26 we show the joint plot of the corresponding real and complex curves. The estimated 



complex (red) curve initially dominates the estimated real (blue) curve (cf. pj'., Fig. 1]). 



1. Close resemblance to generic rank-4 Dyson-index pattern 



It appears now-as a plot (Fig. 27) parallel to that displayed in Fig. 25 indicates-that 
the Dyson-index pattern continues to hold for the range C™^'^"^ G [^,1] in the two-qubit 
generic rank-3 cases, but with the replacement of (2 — 2(7™^"'^^'^)^ by (2 — 2C™"j^"^)^. 



To test the possible applicability of the rank-3 version of the univariance hypothesis (14), 



we estimated the real and complex rank-3 two-qubit Hilbert-Schmidt separability probabili- 



ties using the ESFs displayed in Fig. 26 The values we obtained were 0.208172 and 0.104852, 
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FIG. 27: The rank-3 counterpart to Fig. 25 with the rather well-fitting function (2 — 2C, 



rank—3\ 



replacing (2 - 2C^^^-y. 

respectively (while the correponding conjectured values were perhaps somewhat disappoint- 
ingly different, calling for further analysis, that is ~ 0.235294 and ^ ~ 0.121212). We 
can express these results as one-dimensional integrals over C"^^^'^ G [0, 1] of the product 



of the real function displayed in Fig. 26, and (using C = Cl^'^^ ) the univariate marginal 



max 
rank—3\ 



Hilbert-Schmidt probability distribution (Fig. [28| 

1792^4 (12^2 - 5) 



marg^,^;(C) 



81 



0<C <h 



(20) 



and the integral over C™^ ^ G [0, 1] of the product of the complex function displayed in 



Fig. 26, and the univariate marginal Hilbert-Schmidt probability distribution (Fig. 28) 

™°iC - 1)^C2(C(C(C(16325C - 7693) - 379) + 315) +45) ^ < C < 1 



729 

(155C^ + 1287C^ - 1089^2 + 231) < C < ^ 



(21) 

(These distributions are "marginal", in the sense that they are obtained by integrating the 
HS or Bures measure defined on the three-dimensional simplex of eigenvalues-obtainable 
from the papers of Zyczkowski and Sommers [131 Ej-over two of the three coordinates 
[the third coordinate being Cmax] used to parameterize the simplex.) Also, we have for 



0<C <l, 



marggj^ijjgj,j^jQj^jg (C) 
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(22) 



9209200(7^3 (3 (7133C10 + 236790C8 + 253023C6 - 729980C^ + 497097C2 - 142766) + 46189) 

531441 ' 
(The quaternionic expression for ^ < C < 1 is somewhat more cumbersome in nature to 

present.) To obtain these univariate functions, we have transformed one of the eigenvalues, 

say Ai to C'^2x~^ = ^ -^3 (the jacobian of the transformation being unity) and integrated 

(restricted to the Weyl chamber of ordered eigenvalues) the corresponding (bivariate in this 

case) Hilbert-Schmidt measures (over the eigenvalues) [13, eqs. (4.1), (6.5), (7.8)] over 



A2. Fitting the means and variances of ((20)-(22)), we can obtain beta distribution B{p,q) 



approximations to the real, complex and quaternionic probability distributions using the 
paired sets of parameters {p, g} = ^} ^ {6.62048, 5.73888} , , '^^} ^ 

{10.7838,8.93514} and {SHif , ^^^Hf } - {19.1521,15.3137}, respectively. 
Beta distributions, defined over the unit interval, are a general type of statistical distribution, 
related to the gamma distribution, and have two free parameters. 

The Hilbert-Schmidt (total-separable and nonseparable) probability that a minimally- 
degenerate two-qubit state has maximal concurrence within the range [^, 1] is || = (|)^ ~ 
0.604938 for real states, ff^ ^ 0.667405 for complex states, and ^^^^^J^^^^f^^^^ ^ 0.748123 



for quaternionic states. (For the [smaller] full-rank counterparts see sec. IX D ) In Fig. 29 



we show the counterpart to Fig. 28 based on the Bures (minimal monotone) metric. 

Let us assume (cf. Fig. 27) that the ESF in the real case is proportional over C^^^~^ G 



[|, 1] to (2 — 2CJ^^J:~^)i and in the complex case to the square of this. Then, we have that 
the contributions over this half-domain to the Hilbert-Schmidt real and complex separability 
probabilities, respectively, are 36.5.13.3X ~ 0.380924 (multiplied by a normalization constant 
approximately 0.177365) and 29^39^5.ii^23-2'9^3i ~ 0.327832 (multiplied by a normalization con- 
stant approximately 0.086232). 



F. Rank-five complex qubit-qutrit case 

For the full-rank qubit-qutrit case, the counterpart-although not enjoying all the 
properties-of the two-qubit maximal concurrence formula Q is [24, p. 16] 

CTat"' = max(0, Ai - As - 2 v^), (Ai > A2 > A3 > A4 > A5 > Ag), (23) 

which, obviously (since Ae = 0), degenerates (using the same eigenvalue-ordering) to 

QTr' = Ai - A5. (24) 
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rank- J 



FIG. 28: Univariate marginal probability distributions (20) and (21) over C, 



yrank—3 



of the Hilbert- 



Schmidt measures on the eigenvalues. The quaternionic curve has the highest peak and the real 
curve, the lowest. All three curves are asymmetric about C™^'^~^ = ^ and skewed to the right. 
The real peak (mode) is at 0.534989, the complex at 0.549857, and the quaternionic at 0.558738, 
while the medians are 0.53707, 0.5483 and 0.556787, respectively. Also, the means are ~ 



0.535665, p ^ 0.546875, and |^ 0.555684, respectively. 



27313 
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FIG. 30: Qubit-qutrit eigenvalue-parameterized separability functions as a function of maximal 
concurrence for the generic rank-5 real and complex 6x6 density matrices. 104,200 TF-points 
were used in the real case and 181,600 in the complex case, while, for each such point, the interval 
[0,1] was sampled at 1,000 equally-spaced locations of C'!^2x~^- '^^^ complex (red) curve is initially 
dominant. There are manifest discontinuities (jumps) at 



^rank—5 
^max 



I and 2 in the real case. 



In Fig. 30 -again under the hypothesis (ansatz) that the corresponding eigenvalue- 
parameterized-separability function is a (univariate) function of the maximal concurrence 



expression (24)-we show the analogue of Fig. 26 for the minimally-degenerate (rank-5) 
generic real and complex qubit-qutrit case. (Again, in the complex case we used the SU (6)- 
Euler-angle parameterization of S. Cacciatori P App. A], while we used a yet unpublished 
Euler-angle parameterization of his of S0{6) for the real case.) There are evident jumps 
in the real (blue) curve at C^^^~^ = | and ^. (The still erratic nature of the complex 
[red] curve-we used 1,000 [not 500] equally-spaced points in [0,l]-makes it, at this point of 
sampling, difficult to gauge the applicability of the Dyson indices.) To test the applicability 



of the rank-5 version of the univariance hypothesis (14), we estimated the real and complex 



rank-5 qubit-qutrit Hilbert-Schmidt separability probabilities using the ESFs displayed in 



Fig. [30| The values we obtained were 0.097232 and 0.0226654, respectively, while the corre- 
sponding conjectures were-again somewhat disappointingly different-^ ~ 0.0751174 and 



^ ^ 0.0133445. 
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G. Full-rank real and complex qubit-qutrit cases 



In the two-qubit case, we had evolved a computational strategy in which we used the 
Mathematica command Findlnstance to systematically generate random sets of three or 



four eigenvalues that yielded values of the maximal concurrence or (19)) at equally- 



spaced inteiYaAs Cmax £ [0, l]-and, similarly, in the rank-5 qubit-qutrit case (24). However, 
due to greater complexity in the rank-6 case, this strategy did not prove at all feasible 
for generating random sets of six eigenvalues yielding equally-spaced values of the maximal 



concurrence (23). 

So, we altered our approach, now simply randomly generating density matrices (again 
using the same quasi-Monte Carlo routines [20]) and recording their associated values of 
concurrence. We "binned" these concurrence values into intervals of length ^, and averaged 
the total measures recorded by the number of observations within the individual bins. We 
interpolated these average values to obtain the associated eigenvalue-parameterized separa- 
bility functions (ESFs). We have generated the corresponding curves for both the full-rank 
real and complex generic qubit-qutrit cases, but they are still somewhat crude/rough in 



character. Nevertheless, we plot in Fig. [31] (cf. Figs. 23 and 27) normalized forms of the 
complex (red) curve and the square of the real (blue) curve. They appear to indicate possi- 
ble adherence to the Dyson- index ansatz, since the two curves closely "track" each other, at 
least (as our generally observed pattern in the two-qubit case would suggest) for the higher 
values of Cl^2x~^ ■ seems that this domain of possibly strict Dyson- index behavior may 



be C!^^^~^ G [|,1], while in the full-rank two-qubit case (Fig. 23) it highly convincingly 
appeared to be C™^^~^ G [^,1]- Our level of binning is perhaps too coarse for the detec- 
tion of possible discontinuities in the two curves.) We were interested in seeing how close 
the plotted curves came-under the rank-6 qubit-qutrit version of the univariance hypothesis 



(14)-to yielding the conjectured HS real and complex separability probabilities of ^ and 
[3 sec. X], but the requisite numerical integrations proved quite problematical to per- 



form. In Fig. 32 we plot the two functions in Fig. 31 over the interval [|, 1], along with the 
interweaving curve (| — C^^^~^y. 
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FIG. 31: Joint plot (cf. Fig. 23)~to test Dyson-index ansatz-of numerical estimates of the real 

(The complex [red] curve is lower at C^^J:~^ 



1 / ^rank — Q\ 



and complex 



2 ( f^iVank — Q 



0.) There 



were 62,086,051 20-dimensional TF-points (each point corresponding to a single density matrix) 
used in the real case and 451,373,489 35-dimensional TF-points in the complex case. Each of these 
points was allocated to one of fifty bins in the interval C™„"j=-6 G [0, 1]. 
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IX. SEPARABILITY PROBABILITY DECOMPOSITIONS OVER Cmax RE- 
GIONS 

A. Cmax = domain 

One can-in an apparently natural manner-consider the two-qubit real, complex and 
quaternionic Hilbert-Schmidt separability probabilities to be the sum of three components: 

(1) the Hilbert-Schmidt absolute separability probabilities (corresponding to Cmax = 0); 

(2) the probabilities over the range Cmax £ (0, |]; and (3) the probabilities over the range 
Cmax G [|, !]• (For a contour plot of the three-dimensional body Cmax = 0, see [SI Fig. 2].) 
Now, we have previously been able to compute the absolutely separable components [8], eqs. 
(34), (35)]. These are 

ffo 69 28 - 22 57r , , 

pS±=o = ^ 0.0348338, (25) 

22 



pHb = ^ ^ . — ~ 0.0036582630543035 (26) 

C/maa:— U 016Q5 ^ ' 



,i^5_..e. V'l - ^2^2 - ^3V^+ ^4^2 sec-i(3) 

2I63E 

(the Bures [minimal monotone] metric analogue being considerably smaller, 0.000161792 [8l 
p. 25]) where 

i)i = 956877309536, ^2 = 781862943168, 



and 



where 



V'a = 746624752335, tp^ = 1990999339560, 



Hs , 13 (01 + 02v/2 + 03^271- 04^2 sec-i(3)) , , 

pHA , ^ vv^i — v _ 0.0000398703, 27 



01 = -806338156306739134839776, 02 = 658857590468226345222144, (28) 

03 = 629162653900414735065195, 04 = 1677767077067772626840520. 
(These are "conjecture-free" results, not dependent on any Dyson-index ansatz. In [8^, eqs. 

H S 

(36), (37)] we gave a considerably lengthier, but fully equivalent, expression for Pc^^^=o- One 
might seek to find explanations for the large integers displayed above in terms of gamma 
functions. The computational challenges to computing analogous absolute separability re- 
sults for the qubit-qutrit states appear to be highly formidable.) 
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Further, accepting the strongly-supported Dyson- index ansatz (Fig. 25 and (p^) that 
o^^\C^a.{^,...X,)) = and a(2)(C„,,.(Ai...A4)) = for ^na. G [^,1], 

we can now add to the absolute separability probabilities {Cmax = 0) listed immediately 
above, the conjectured probability contributions 



P, 



(3162214 - 738885v^tan~^ (v^)) 



Cmax £ [ 2 ;!] 



)12 



5 ■ 7 - 17- 19 



0.02559647778, 



(29) 



and 



7 f 148453588142 - 79729806357^2 cot'^ 
Further, using the Dyson-index ansatz with /? = 4, we obtain 



0.01029059519. 



(30) 



266 . 38 . 11 . 29 . 31 



K ^ 0.165191k, 



(31) 



where 



Ci = 174916374035295022487516506, C2 = 42964561240209557008032951, 



(32) 



and K is the /3 = 4 unknown and yet-unconjectured analogue of the presumed real and 



complex constants -j^ and j^. (In computing (29)-(32), we found a joint transformation of 



the form ai 



^ and a2 = ^/Xq^h to be helpful.) 



1. Corollaries to the "twofold" SBZ-Theorem 



Since the probability is zero that a generic minimally-degenerate two-qubit state is abso- 
lutely separable (that is, P"^"-"^^-^ = P^""^g^Q i])-as can be immediately deduced from (24)-we 



have simple corollaries to the twofold-theorem of Szarek, Bengtsson and Zyczkowski [36] of 



the form 



-yrank—A 

C,„a^e[o,i] 



-prank— 

c„a^e[o,i] 



2 - 



^rank—4 



P. 



rank— 4 

Cmax —0 



prank— 3 

C„a:.e[0,l] 



(33) 



prank— 3 ' 

where the P's are Hilbert-Schmidt separability probabilities, for the real, complex or quater 
nionic two-qubit states. 
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FIG. 33: Previously-generated |17| Fig. 1] estimates of the two-qubit real and complex eigenvalue- 
parameterized separability functions over the domain Cmax £ (0, The complex (red) curve is 
initially higher-valued and the real (blue) curve, close to linear. 

C. Cmax £ (0, 2] 

So, the most conspicuous missing parts in the Hilbert-Schmidt separability probability 
"puzzle" appear to us to be formulas for P^^''^"' 1,, and P^^'""^* 1,. Of course, we 

can subtract the sums of the other two parts (P"^'-^^' + Pf'^''^"' 1 p^^'"^^'-'^ + p^^-™p^- 
and P^'^''"^n + Pif'^''"'"ri n) from our overall coniectures of t^.^^ and J„^iff„^„tt to obtain 
"induced" conjectures about these third components. 

Since it now appears crucial to, additionally, model the eigenvalue-parameterized sepa- 
rability functions over the domain Cmax & (0, we present in Fig. 
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for the convenience 



of the interested reader, the previously-generated [T71 Fig. 1] estimates of these functions. 
The real (blue) curve is close to linear (^ 1 — 1.75Cmax)- Also, we have noted that the 
complex (red) curve is quite well- fitted by cos^^ (Cmax) and cos^ {2Cmax)- However, we ap- 
pear here to lack a strictly similar Dyson-index ansatz to serve as a guide in constructing 
these two functions (cf. [2|). Also, there were indications given in pTl Figs. 3-5] that these 
two functions have multiple (matching) points of discontinuity in Cmax ^ (0,^]- (These 
were Cmax ~ 0.204, 0.294, 0.34.) The lack of a Dyson-index pattern strictly similar to that 
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FIG. 34: Estimated ratio (for the full-rank two-qubit case) of 



2 ( j^rank — 4:\ 



along with the approx- 



imating lines 1 -|- 2Cmax and 2 the vertical lines C„ 



^ and 1 
100 2- 



found apparently for Cmax ^ l^A] to exploit for Cmax ^ (0, ^] is immediately apparent from 
A flat line over at least some subdomain of C^ax G (0, ^] would indicate such a 



Fig 
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Dyson-index pattern. Clearly, no such flatness appears there. However, it now seems that 
there is a pattern of the approximate form 

2 ■ • 1, 



2 / /^rank—4\ 
y^max J 



rank— 4'' 
max 



rank— 4 



e (0, 



(34) 



(We have already noted that cr^(C™"j'~'^) ^ 1 — l.TSC^^j""^ in this half-domain.) The highly 
interesting nature of Fig. [34] led us to similarly re-examine the minimally degenerate rank-3 

serves as an 



two-qubit scenarios (Fig. 26). Thus, we obtained Fig. 35 Now 1 + 3C. 



excellent linear approximation, and we have a relation analogous to (|34|) 

I t '2//-irank — 3\ 

yrank-4:\\2 ^1^^ t (L^max ) 



yrank—3' 
''max 



2 T 

^rank—3 ' - ^ 



e (0, 



(35) 



i^KC:::t')y and also 



-I .n/^Tank — 3 



^HC:::t')f over the 



2 /^ranfc — 4\ 

If we plot " ^^.T/nfe-i vs. 
half-domain, the two curves within each set are essentially indistinguishable. We have 
investigated analogous plots of the same form for the minimally- degenerate qubit-qutrit 



(Fig. 30) and full-rank (Fig. 31) cases over Cmax £ [0, \]. They are much rougher in nature. 



due to our limited sampling, but still indicate initial monotonically-increasing (non-flat) 
behavior over Cmnr- 
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4 




^2 (fj'rank — 3\ 

FIG. 35: Estimated ratio (for the minimally-degenerate rank-three two-qubit case) of ^ Znk~3-.\2 



along with the approximating lines 1 -|- ^Cmax and M and the vertical lines C, 



and 



D. Total probabilities over Crnax regions 

Let us also point out that the Hilbert- Schmidt probability (P) that a generic two-qubit 
real state (separable or entangled) lies in the domain Cmax ^ [^5 1] is 

75962 - 25515v^tan-i (v^) 



pH Sreal 

with the complex counterpart being 
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0.187584, 



(36) 



p: 



174957361466 - 124912178055v^cot-i (^^ 



and the quaternionic analogue being 



231 . 35 



0.241961, 



(37) 



Cmax ^ [ 2 ' 



7i-72V^cot-i {£) 

263 . 310 



0.323053, 



where 



7i = 217894901318574565900294, 72 = 107614737772623370233945. 



(38) 



(39) 



(The comparable total probabilities for the minimally-degenerate two-qubit states have been 
given in sec. VIII E 1 ) Since for the absolutely separable states {Cmax = 0), the two proba- 



bilities Pcmax=o -Pc™a:r=o Sltq equlvaleut, we can immediately determine (by subtracting 
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from 1 our known HS probabilities Pc^^^=o and Pcmax&[- i]) complementary probabili- 
ties, Pr prn Numerically, these are P/f^''^"' i, ^ 0.777582, P^^^°*, f« 0.754381 and 
P/^^^""* 1, ~ 0.676907. 



X. CONCLUDING REMARKS 

Our analyses of two-qubit dmg'ona/- entry-parameterized separability functions (DESFs) 
[3 [311 SS] cind ez^'enwa/we-parameterized separability functions (ESFs) [U [17] completely 
share a common goal: the determination of two-qubit separability volumes and probabilities 
(in terms of various metrics). As pieces of these formidable objectives begin to be assembled, 
we can pose a further challenge-to find transformations between the two different sets of 
coordinates used-that is, (1) the diagonal entries and (2) the eigenvalues of 4 x 4 density 
matrices-that will map one set of separability functions into the other. The Schur-Horn 
Theorem, which asserts that the decreasingly-ordered vector of eigenvalues of an Hermitian 
matrix majorizes the decreasingly-ordered vector of its diagonal entries |l6l chap. 4] (cf. 
[m [IE]), would appear to be of possible relevance in this regard, particularly since the 



maximal concurrence Cmax over spectral orbits (12) is expressed in terms of the ordered 
eigenvalues. 

In terms of the diagonal entries {Di, D2, -D3, D4 = 1 — Di — D2 — -D3) of 4 x 4 two-qubit 
density matrices, we can express the conjectured Hilbert-Schmidt separability probability 
[7], Jj, of generic complex states in the form 



12,108,096,000 




(P'lP'aP'aA) (3 - uyudDidD2dD3, (40) 
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where u = and the integration extends over the unit simplex, but with the restriction 

u <1. (We note that 12, 108, 096, 000 = 2^ ■ 3^ ■ 5^ ■ 7^ ■ 11 ■ 13. Let us also observe that 
the variable log u conveniently ranges over the entire real axis and is symmetric about the 
origin.) 

Additionally, in terms of the eigenvalues (Ai, A2, A3, 1 — Ai — A2 — A3) of 4 x 4 two-qubit 



density matrices, we can express this same separability probability as (cf. (14)) 
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2,201,472,000 / / / or(2)(C„,,(Ai . . . A4))nf (A, - A.O^rfAA^/Aa, (41) 



and the integration extends over that part (Weyl chamber [16]) of the unit simplex for which 
Ai > A2 > A3 > A4. (We note that, interestingly, in light of the just previous factorization, 
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2, 201, 472, 000 = 2^° ■ 3=^ ■ 5^ ■ 7^ ■ 13.) Here a(2)(C^„^(Ai . . . A4)) is the (two-qubit complex 
[/3 = 2]) eigenvalue-parameterized separability function that we have previously sought to 



determine [T71 Fig. 1], and was found to be very well-fitted by 



{2 SCjTiaa;)^ 

15 



for Cmax e [^,1] 



(Fig. 25 and (18)). (It is possible to reexpress these two last integrals so that both are taken 



over the same complete 3-dimensional unit simplex.) Further still, our generic complex 
two-qubit Bures separability probability conjecture ^ [H Table VI] takes the form 



1680(\/2 - 1) 



TT 



7r= 




a^'KCmaA\i---X4))^, (A. -A,) 



n-< /\' , d\^d\2dXs. (42) 



71680 y J J VA1A2A3A4 "'^^ + 
It is abundantly clear: (a) that this (piecewise continuous) function (J^'^\Cmaxi^i ■ ■ ■ A4)) 
has a jump discontinuity at Cmax = | (as well as does its real counterpart 
c"''^'' (C'maa; (Ai ... A4)); and (b) that in the diagonal-entry-parameterized scenario, the value 
u = = 1 is a locus of special symmetry. In this regard, we might speculate that if one 
can find a coordinate transformation between the two separability probability expressions 
((40) and (41)), then those values of the Aj's for which Cmax = ^ wiU be mapped to those 
values of the Dj's for which u = 1. 

Through the use of the jacobian transformation of the diagonal entry (say) to u [I5l 
eq. (11)], and subsequent integration over Di and D2, it is possible to explicitly reduce the 



computation of the trivariate integral (40) to that of a univariate integral in z/. In Fig. 36 



we show-based on numerical calculations-the univariate marginal probability distributions 
of the Hilbert-Schmidt measure over the real and complex two-qubit states in terms of 
^max~'^ (cf. (20)-(22)). Similarly to their rank-3 counterparts( Fig. 28), these curves have 



differently-positioned peaks, and are not symmetric, but skewed to the right. (We take the 
range of C™^^~^ to be [—^,1], to accord with actual values, rather than the conventional 



[0, 1] [cf. ([12|].) In Fig. 37, we show the Bures-metric counterpart, although we encounter 
some "glitch" in displaying the real curve here. 



The counterparts to the formulas (40) and (41), in light of our conjecture [7] that the 
Hilbert-Schmidt separability probability of the generic real two-qubit states is ^, are (the 
domains of integration being the same) 
1,209,600 



8 
17 



17 



J J J{D,D2DsD^f/\3-u)y^dD,dD2dDs, (43) 



and 



8 15,482,880 



17 



17 




a^'\Cmax{Xl ■ ■ ■ A4))Hf<,(A, - X,)dX,dX2dX 



(44) 
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-a.4 -CiJ2 02 0.-^ 0.6 Ci.S 3 JO 

FIG. 36: Full-rank two-qubit univariate marginal probability distributions over C'l^^^~^ € [— |, 1] 
of the three-dimensional Hilbert-Schmidt measures on the eigenvalues. The quaternionic curve has 
the highest peak, and the real curve, the lowest. Numerical methods were used. 
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FIG. 37: Counterpart to Fig. 36 based on the Bures metric. The quaternionic curve has the highest 
peak, followed by the complex curve, with the (incompletely constructed) real curve, apparently 
the lowest. Numerical methods were used. 
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(Here, 1, 209, 600 = 2^ ■ 3^ ■ 5^ ■ 7 and 15, 482, 880 = 2^^ ■ 3^ ■ 5 ■ 7. It appears [Fig. [27| and 



hsh that possibly a^'KCma.{Xi . . . A4)) = for C^a. G !]•) 



Let us point out the possible relevance of the concept of the Thouless energy [IHl p. 
734] in the modeling of the threshold or crossover effect we have numerically observed for 
eigenvalue-parameterized separability functions in both the full generic real and complex 
two-qubit and qubit-qutrit cases. There, the Dyson indices (/? = 1, 2) of random matrix 
theory only seemed to apply above a certain value of the maximal concurrence Cmax (that is, 
I in the two-qubit case, and possibly | in the qubit-qutrit instance). (It remains to formally 
reconcile these observations with the ones that, in terms of diagonal-entry-parameterized 
separability functions, Dyson-index behavior appear to be strictly followed [7].) 
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